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Abstract: We work out the optimization problem, initiated by K. Soundararajan, for 
the choice of the underlying polynomial P used in the construction of the weight func- 
tion in the Goldston-Pintz-Yildirim method for finding small gaps between primes. 
First we reformulate to a maximization problem on L?[0, 1] for a self-adjoint oper- 
ator T, the norm of which is then the maximal eigenvalue of T. To find eigenfunc- 
tions and eigenvalues, we derive a differential equation which can be explicitly solved. 
The aimed maximal value is S(k) = 4/(k + ck!?), achieved by the k — 1“ integral of 
xl K/27 (a, Vx), where a, ~ ck!? is the first positive root of the k — 2" Bessel func- 
tion J, ». As this naturally gives rise to a number of technical problems in the appli- 
cation of the GPY method, we also construct a polynomial P which is a simpler func- 
tion yet it furnishes an approximately optimal extremal quantity, 4/(k + Ck!) with 
some other constant C. In the forthcoming paper of J. Pintz [8] it is indeed shown how 
this quasi-optimal choice of the polynomial in the weight finally can exploit the GPY 
method to its theoretical limits. 


Keywords: Prime distribution, gaps between primes, Goldston-Pintz-Yildirim method, 
self-adjoint Fredholm type operator, Bessel differential equation, Bessel functions of 
the first and second kind 


Classification: Primary 11NO5, 47475. Secondary 47453, 49]05, 49K15, 49N10 


Balint Farkas: Bergische Universitat Wuppertal, Faculty of Mathematics and Natural Science, Gaufss- 
trafe 20, 42119 Wuppertal, Germany, E-mail: farkas@uni-wuppertal.de 

János Pintz: Hungarian Academy of Sciences, Alfréd Rényi Institute of Mathematics, 1364 Budapest, 
P.O. Box 127, Hungary, E-mail: pintz@renyi.mta.hu, revesz.szilard@renyi.mta.hu 

Szilárd Révész: Kuwait University, Department of Mathematics, P.O. Box 5969 Safat, 13060 Kuwait, 
E-mail: szilard@sci.kuniv.edu.kw 


1 Introduction 


1.1 The extremal problem as given by Soundararajan 


In his work [11] Soundararajan presents and analyzes the proof of Goldston-Pintz- 
Yildirim yielding small gaps between primes. Among others he raises and answers one 
ofthe most important problems of the field: Is it possible to modify the weight function 


Bereitgestellt von | De Gruyter / TCS 
Angemeldet 
Heruntergeladen am | 16.10.19 13:23 


76 = Balint Farkas, Janos Pintz, and Szilárd Révész 


a(n) in such a way that the method would lead to infinitely many bounded gaps be- 
tween consecutive primes. If we consider the weight functions in full generality, that 
is all functions a(n), then this leads essentially to a tautology. For example, defining 
a(n) = 1if both n and n + 2 are primes, and otherwise setting a(n) = 0, the summatory 
function of a(n) describes the number of twin primes up to x. Thus we cannot hope for 
an asymptotic evaluation of the summatory function. We briefly describe the feasible 
choices of the weight function a. Let us take an admissible k-tuple H = {h,,...,h,} 
meaning that there is no prime p with the property that the elements h; of H cover all 
residue classes mod p. Let P(n) = II. * hj) and let us define A, = ud P(S) 
with a nice function P, for example a polynomial, with the additional property A, = 1, 
which is equivalent to P(1) = 1. Afterwards we reduce our choice of a(n) to those of 
type a(n) = Y 4n ap, ( ^4 and try to evaluate the summatory function of a(n) and that 
of a(n)x(n + h), where h is an arbitrary number with h < log n and y is the characteris- 
tic function of the primes. (In case of bounded gaps between primes it is sufficient to 
consider the case when h = h; i = 1,2,...,k.) 

Soundararajan explains how the optimal weight function a(n), hence A4, should 
be chosen to obtain best result: see formula (8) in [11]. In order to get this optimum, he 
also explains the choice A, := wa) P(E) where P is some suitably nice function, 
like a polynomial or at least a sufficiently many times (at least k times) differentiable, 
smooth function on [0, 1] (or at least on [0, 1)), vanishing at least in the order k at 0, and 
satisfying the normalization P(1) = 1. Then, according to the analysis by Soundarara- 
jan, the optimal choice for a(n) and A, is equivalent to looking for the maximal possi- 


ble value of (12) of [11], i.e., to determining 


k-2 


1 1 k 
"E x (k-1) E 2 x (k) H 2 
S(k) :- (| os (PPa - x)) ax) /( | a (PPa - x)) dx), (1 


0 


where the set of functions P, to be taken into account in the supremum, can be the 
set of certain polynomials as before, or more generally a family of functions subject to 
some conditions. 

Soundararajan [11] shows that the question whether we are able to find in this way 
infinitely many bounded gaps between primes is equivalent to the problem whether 
there exists any natural number k with S(k) > 4/k. Then he mentions that the opposite 
inequality S(k) « 4/k holds for all k and therefore the method cannot yield infinitely 
many bounded prime gaps. (In an earlier unpublished note [10] he gives the short 
proof of this fact; we will reproduce this in Subsection 2.3. His considerations also lead 
easily to the stronger inequality S(k) < 4/(k + c log k), cf. Subsection 2.3). Although his 
work answered negatively the above mentioned central problem, it gave some hints 
but did not answer the question: What is the best weight function that can be chosen, 
and what size of gaps are implied by it? In their work [5] Goldston, Pintz and Yildirim 
showed that if one takes P(x) = x**’, where k and £ are allowed to tend to infinity with 
the size N of the primes considered, then with several essential modifications of the 
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original method one can reach infinitely many prime gaps of size essentially log p. 
(To have an idea of the difficulties it is enough to mention that the rather condensed 
proof ofthe result needs about 40 additional pages beyond the original one, presented 
with many details and explanations in [4]. However, a shortened, simplified and more 
condensed version [6] needs only 5 pages). In this case £ = c Vk and the value of the 
fraction (1) is 4/(k + c' Vk) for the given choice of P(x) = xE, Beyond the mentioned 
important fact that k and £ are unbounded in [5], the scheme of the proof is similar 
but not the same as in the simplified version of Soundararajan [11]. However, a careful 
analysis suggests that in order to find the limits of the method it is necessary (but as 
discussed a little later, not necessarily sufficient) to find the size of S(k) as k tends to 
infinity together with the function P which yields a maximum (if it exists) in the supre- 
mum, or at least a function P which yields a value “enough close” to the supremum. 


1.2 Conditions and normalizations 


Before proceeding, let us discuss right here the issue of conditions and normalizations 
in the formulation of this maximization problem. First, it is clear that P“ remains un- 
changed, if we add any constant to P^-?, Thus the extremal problem becomes un- 
bounded under addition of a free constant, hence at least some conditions must cer- 
tainly control this divergence. 

In the number theory construction of Goldston—Pintz—Yildinm, (by now gener- 
ally abbreviated as the “GPY method") the natural restriction is that P must be a poly- 
nomial divisible by x* — or, if we try to generalize the method, then a k-times con- 
tinuously differentiable function with pO vanishing at 0 for j = 0,..., k — 1. That is 
P(x) = PEO xk + o(x**!), The reason for that is the fact that the whole idea hinges 
upon the use of the generalized Mobius inversion, more precisely of the A ; function, 
which must be zero for numbers having at least k prime factors — always satisfied by 
the numbers represented by the product form (n+h,) ---(n+h,) in the construction. So 
for any meaningful weight function we need to use weights not containing any smaller 
power x! than x“. In other words, we should assume here P having a zero of order k, 
i.e., P(0) = P'(0) =--- = p*-P(9) = 0, while P (0) can be arbitrary. 

The analysis of Soundararajan exposed the question, whether a linear combina- 
tion of monomials, i.e., a polynomial, or perhaps some more sophisticated choice of a 
weight function, may perhaps improve even upon this. We can say that the theoretical 
limit of the GPY method is the result, obtainable in principle by a choice of the weight 
function P maximizing the extremal quantity (1). Yet it is to be noted that the techni- 
calities of GPY are far more substantial than to simply “substituting any P" in it would 
automatically lead to a result — it is not even that clear, what result would follow from 
a given weight function. Therefore, to test the limits of the GPY method, we should 
break our approach into two parts. First, we look for the optimization of the weight P, 
in the sense of (1), and second, we extend the GPY method using that weight function. 
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This paper is concerned with this first question, and the second part of this program 
is carried out in [8]. 

The aim of the present analysis is to settle the issue of optimization in problem (1). 
We find the optimal order, and the maximizer of the problem (1), furthermore, as this 
maximum can be achieved by a relatively sophisticated choice of the weight function 
P — actually a transformed Bessel function — we also construct a polynomial weight 
which is approximately optimal in (1). 

Part of these results were reached by J. B. Conrey and his colleagues at the Amer- 
ican Institute of Mathematics already in 2005. Using a calculus of variation argument 
they found the Bessel function J,_, and made some calculations for concrete values of 
k (without analyzing the case k — oo). The fact that the Bessel functions may perform 
better than polynomials in the GPY method is also briefly noted in the book of J. B. 
Friedlander and H. Iwaniec [3] without going into details. 


1.3 Structure of the paper 


In this paper we proceed along the following course. 

Interpreting the problem in the widest possible function class which makes sense 
(i.e., when at least the occurring integrals exist finitely), in Section 3 we make several 
further reformulations until we arrive at a maximization problem in the Hilbert space 
L?[0, 1]. Exploiting the rich structure of Hilbert spaces, and the particular properties 
of the reformulation as a certain quadratic form with a Fredholm-type operator, we 
derive existence of maximizing functions in this wide function class. Then we also 
exploit the concrete form of the kernel in our Fredholm-type operator and compute 
that the maximizers, or, more generally, eigenfunctions, are necessarily smooth. Fur- 
thermore, in Section 4 we find that these eigenfunctions satisfy certain differential 
equations. The solutions are then found to be transformed variants of certain Bessel 
functions. Also it turns out that the solutions are analytic, and they yield a function 
value in the extremal problem directly related to the choice of a parameter, which, 
due to the initial value restriction P*-?(0) = 0, must be a zero of the arising Bessel 
function J,_,. Finally these combine to the full description of the maximal value S(k) 
together with the precise form of the extremal function. From the well-known asymp- 
totic formula for the first zero of the Bessel function J,,, when m — oo, we derive that 
S(k) is precisely asymptotic to — with a concretely known constant c = 3.7115.... 

Unfortunately, in spite of analyticity and power series expansion, the found ex- 
tremal function is too complicated to be used in the number theory method of GPY. 
Basically, we need restrictions on the degree and the coefficient size for the powers 
appearing in the weight function P to make the complicated method work in a tech- 
nically feasible way. As discussed in Section 5, not even calculations using the power 
series expansion of Bessel functions lead to feasible expressions. Therefore, finally we 
look for quasi-optimal polynomials, which still achieve close to extremal values. The 
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result of the last section is the concrete construction of a polynomial P pausa the 
needed technical requirements and still achieving in (1) a ratio of the order = au z with 
some other constant C. That suffices in the method of GPY, because the value of the 
constant C does not increase the order, only the arising constants, in the final result. 

Settling the issue of the search for optimal and quasi-optimal weights, the door 
opens up for revisiting the method of GPY and not only improving upon all the known 
results, but also pushing the available techniques to the theoretical limits of that 
method. This closely connected work is carried out in the paper [8]. 


2 Reformulations and the finiteness of S(k) 


2.1 Reformulations 


The normalization P(1) = 1 is rather inconvenient because the next reformulation 
(still following Soundararajan) is to put Q(y) :- pi y), a completely logical step 
in view of the fact that no values of P, P’, etc. P? occur in the actual optimization 
problem (1) and that the still occurring P7? and P? can be nicely expressed as Q and 
Q'. So in line with the restriction that P vanishes at least to the order k at 0, following 
Soundararajan we write 


P(x) = i 
Uh 


1 Xk-2 


] pe P dx, ...dx, = | | vee | P (ds, , e dx, 
0 0 0 


1 Xk-2 Xi Xk-1 


| Q(x, ,)dxy ,...dx, = ite | Q'(xi)dx, ... dx. 
0 0 


C —— k Ow X 


Therefore, P(0) = -.. = P2 (0) = P*®YD(o) = 0 transforms to the simpler requirement 
that Q(0) = 0, while the corresponding P is obtained by the above integrals directly. 
Let us record one more thing here: The condition that P(1) = 1, expressed in terms of 
Q, is a linear restriction, as I(Q) := (P(1) =) IN j n jaa Q(x, dx, , ...dx, is justa 
linear functional on the function Q. To express it in a more condensed, closed form, 
we may apply Fubini's theorem to get a representation in the form of the well-known 
Liouville integral 


(x —1)*? 
P(x) = jos dor ^ 
1 1 
. 1- De -2 k-2 
P(1) = | Qr ort [ou -»g- z“ 


Note the similarity to the numerator of the quotient in (1). It is thus immediate by 
the Cauchy-Schwarz inequality that P(1) is a finite, convergent integral whenever the 
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Lebesgue integral if Q’(1- y) dy exists. That is, no special requirement is needed 
to this effect once we guarantee that the numerator and denominator in (1) are well 
defined. 
In all, we were to look for maximum in the family 
x Xi Xk-2 
piu [ren - | | E | pU dx, ...dx,, P(ü)-1 P*%(0) = of (2) 
00 0 
but following Soundararajan we changed the setup to 


l y? 
Q; := fo : Jaa = NEY =1, Q(0- o}. G) 
0 


where now Q can be any (say, continuously differentiable) function satisfying the re- 
quirements. This also means that we want the occurring functions to belong to a suit- 
able function class, to be specified later. The quantity we seek to maximize is then 
expressed as 


T ous 
n ny" !Q*(1- y)dy 
1 
h "Q^à - y)dy 
which again is a fraction of two expressions, both quadratic homogeneous in Q. There- 


fore, the ratio will be the same for cQ with any c + 0 and the original question can thus 
be rewritten as looking for the supremum of these quantities among functions in 


(n:2 k- 1), (4) 


1 
Q* := lo : [ea - y)dy #0, Q(0) = o}. 
0 
Continuity of Q' is not indispensable, but of course the ratio must be a ratio of finite 
quantities, with a non-zero and finite denominator, hence we need still to restrict con- 
siderations to functions Q' £ 0 or, in general allowing discontinuous functions, Q' not 
zero almost everywhere and also satisfying th x"Q"(1- x)dx < oo. 

We will see in a moment - see the proof of the forthcoming Proposition 1 - that 
this latter condition also implies that even fax” Va — x)dx < oo, as needed. Fur- 
thermore, together with the restriction that Q(0) = 0, we see that Q is constant if only 
Q = 0, so we need to exclude only this obviously singular case. Otherwise also the 
numerator remans finite, i.e. the ratio (4) exists finitely, whence S(k) exists at least as 
a supremum of certain finite, positive quantities. 

Let us observe that the condition that P(1) # 0, is a linear condition, equivalently 
stated in the form that the linear functional! Q — IN Q(1 - y)y" ! dy on the function 


1 Linearity is clear, once the integral is defined finitely. Then again, finiteness of IR ya - y)dy, 
appearing in the numerator of the extremal quantity, ensures by means of the Cauchy-Schwarz in- 
equality, finiteness of the functional values, too. So for the rest of the argument to be valid, it suffices 
to check finiteness of the numerator of (4). 
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space of admissible functions Q should not vanish. In other words, the subset which 
falls out of consideration for not meeting this condition is the kernel subspace of the 
linear functional, which is of codimension one - in view of the fact that the functional 
itself is not identically zero, obvious from looking at functions Q with Qļo,ı) > 0 cer- 
tainly yielding positive functional values - so a hyperplane 2€ of our linear function 
space X (whatever choice of the function space and respective norm we make later 
on). 

Therefore, dropping the restriction that P(1) = Q(0) + 0 means only that q € KH 
is dropped. In the following we will find the supremum on X, and actually will show 
that here the supremum is finite, attained at certain maximizers. 

The only issue which may bother us a little, is if the actual maximizers will belong 
to the singular hyperplane «K, or stay in X V H. That we should check at the end. But 
maximizers q € X will actually be positive functions, so the value of the functional 
I(q) = P(1) will be necessarily positive for these, and maximum over X or X \ H will 
thus be seen to be the same. We will leave it to the reader to check this and from now 
on pass on to the class X. 


2.2 The choice of the function class of the extremal problem (1) 


In view of the above, let us fix the function class, where the extremal problem (1), 
initiated by Soundararajan, will be considered, as follows. Write q(x) = Q(1- x) as 
before. Then the whole problem becomes 


1 1 
max (k - y Jo) / ( J^ do) under condition q(1) 2 0, (5) 


0 0 


understood in an appropriate function class X, like, e.g., C'[0, 1]. 
Partial integration in the numerator and q(1) = 0 yields now the reformulation 


1 1 
S(k) = sup ( -2 | sigaar) / ( | awa), (6) 
qEX,q(1)=0,4+0 


0 0 


Certainly we want the denominator to be finite, so we assume that our function class 
is chosen in such a way that for any q € X this weighted square integral of q’ con- 
verges. This implies the convergence of the numerator (as we'll see soon) and conse- 
quently that also the positive, non-degenerate linear functional I(q) := (k - 2)!P(1) = 
I^ q(t)t*?dt is well defined, finite. So now we fix the largest function space we may 
deal with as 


1 1 
X := t :(0,1] >R : q(x) =- | q (t)dt, [44^ cod < Jl ; (7) 


x 0 
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The definition above is understood to mean that any q € X is an absolutely continuous 
function on each compact subinterval of (0, 1], whence q € EL (0; 1] and q(x) exists 
as a Lebesgue integral of q’, and in view of the second condition, q' is also square-in- 


tegrable on [0, 1] with respect to the weight x*. 


2.3 An estimation of the extremal value 


Before proceeding let us stop for a little further analysis, establishing boundedness of 
S(k), because this will be needed in what follows. 

Soundararajan [11] remarks that *the unfortunate inequality" S(k) « 4/k holds. 
This is not completely obvious, but in fact the situation is even worse, namely, S(k) « 
4/(k + clog k). This was essentially proved (without an explicit calculation of c) in the 
mentioned unpublished note of Soundararajan [10]. Together with the mentioned ex- 
ample P(x) = x***, £ = c vk, this shows that the value of S(k) is between 4/(k +c’ log k) 
and 4/(k + c” Vk). 


Proposition 1 (Soundararajan). The extremal problem (1) is bounded by 4/k. Moreover, 


we have S(k) < Funes for allk > 4. 


Proof. Let us fix, as in (4), the value n := k - 1. We are to show that whenever the 
denominator of (4) exists finitely, but is non-zero (i.e. when P? = Q' + 0), then also 
the numerator (with the condition that Q(0) = 0, i.e. p-U(g) - 0) exists finitely and, 
moreover, the ratio admits the stated bounds. 

Let us write now q(x) := Q(1 - x), assume that Q, hence also q, are absolutely 
continuous, and consider the resulting relations q (x) = -Q'ü- x), q(1) = Q(0) = 0. 
These imply by absolute continuity that q(x) = q(x) - q(1) = - f. q (t)dt = rh Q'ü - 
y)dy = [-Q(1 -yl = Q(1-x) - so we can as well start with the conditions that p(x) := 
q (x) is measurable and finite a.e., admits the weighted bound L := |: x"q" (x)dx « co 
(coming from the requirement that the denominator is finite), and also that q’ does not 
vanish a.e. (for the denominator being positive). First let us check that then defining q 
from the given p := q' as q(x) := - f. q (t)dt works, results in an absolutely continuous 
function, and with this function the numerator stays finite, bounded in terms of L. 

Indeed, f^ Iq (Dldt < yf} dt f tla (Pat < VE fo mg pac = ax 
(valid for all n > 2) gives not only that q(x), as a Lebesgue integral, exists for all x, but 

L l-n 


also the estimate |q(x) P< =x "on (0, 1]. It follows that q(x) is absolutely continuous 


with derivative p = q' a.e. Moreover, f; nx" g^ (x)dx < —_L, too, hence the numerator 
in (4) is also finite and the quotient cannot exceed = « 2 (k > 3). 

Fork = 2, i.e. n = 1, there is only a little difference in the calculation, as then we 
obtain |q(x)| < VITlogxTand f; nx" q'Gdx = f, q?(x)dx < L f C-logx)dx = L < 2L, 


extending the above bound to all k > 2, too. 
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In the following we compute the stated sharper bound, too. For any m » n - 1, 
I(m) := IN x"q'(x)dx < I x" g'(x)dx < 2L < oo, as by condition we consider the 
class of functions satisfying 0 < L < oo (with L := jix ng? (x)dx). 

Partial integration (using also q(1) = Q(0) = 0) and Cauchy-Schwarz inequality 
yield 


Pu 


I(m) = 2q' (x)q(x)dx < 


2 
i VI(2m +2 — n)L. 


| 
um 
= 


So starting from m := mọ := k — 2 and continuing by induction with m; := k + 2-3 
(j = 0, 1,..., N), we arrive at 


—(N+1) 


f, xq (x)dx dm) | Iu 2 M (ka 25 - 3)? 
fxCuhodx E R22 [ue 


Since 0 < I(v) is decreasing with v, I(k + 2N*! — 3) converges with N so that we can 
pass to the limit N — co, and then even take supremum with respect to q, obtaining 


p 4 4 e PER EA 
sts e-nf1( 2) “Tie &+2-2 G+) ; 
ju 


27 
(k+ -2 k-2y 


Observe that for all j > 1 every single k«2/ -2 > kin the denominator of the last but one 
expression, hence S(k) < 4/k follows immediately. We can even sharpen this estimate 
further. Let us denote the last product by D := D(k) and define £ := [log, (k - 2)]. Then, 
by using log(1 4 x) » x — né (for x » 0) we infer 


oo log(1+ 4) £l ji 1 y 2 
og Dw) - -$ EC Yves) 


Therefore, as e'* < — (for x > 0), we obtain 


"P exp ( Gi 4 4 4 

a d A ee 
k-2 k-2 Kee) k-£-3 = ee log k= 5 
since £ > log, (k - 2) - 1 2 log, k - 2 for all k 2 4. m 


A further elementary observation is that for the Cauchy-Schwarz inequality to 
be precise, we should have x"*!-"/?g = cx"?4$' in all the above applications of the 
Cauchy-Schwarz estimate (i.e., for all occurring values of m). This cannot hold for 
whatever choice of q for all m simultaneously. To have an about optimal estimate we 
may strive for having the Cauchy-Schwarz estimate sharp at the very first application, 
when m = k - 2 andn = k - 1, so q = cq’ follows, and then q(x) = e**. But even that 
is not a valid choice in our problem: q(1) = 0 prevents us taking q as an exponential 
function as it can never be zero. In any case, the estimate of S(k) above cannot be 
sharp. 
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3 Existence and smoothness of maximizers in the 
extremal problem 


3.1 Existence of maximizing functions in the extremal problem 


In this paper the role of k is fixed. Furthermore, it will be convenient for us to avoid 
repetitious use of k - 2 and k - 1, so throughout the rest of the paper except for the last 
section, Section 6, we will fix the notations for two further integer parameters. So we 
define 

m:-k-2, n:-k-1. (8) 


As explained above, we can discuss the optimization problem in the function 
space 
1 1 
y: fo € L1. (0; 1] : [Pora < e where q(x) =- | pdt. (9) 
0 x 


Multiplying the occurring functions by t”, we can even consider the space of func- 


tions g(t) := pne ?^2g (£t, which then will be square-integrable on [0, 1], so that 
ọ € I?[0, 1]. 

Next let us establish, how the functional to be maximized looks like over these 
spaces. On X, on Y and finally on I? [0, 1] we must consider the respective equivalent 
expressions 


-2 b x*!q'(x)q(x)dx F 2 j x” p(x) (f. p(t)dt) dx 
h xg" (x)dx 7 [a p?(x)x"dx 
2 fo ecox"? (f. pe"? at) dx 2f, fo PPO Rox Go Dx" ^E" dtd 
H g? (x)dx k g? (x)dx 


: j. IN qG)g OK (x, t)dtdx with K(x á u pu t) » (10) 
= j; e! G)dx "57 \ max(x, t) , 


the last step being a technical one to bring the kernel K to a symmetric form. So finally 
we find that 
fo fo POPOK, Ddtdx 


S(k)2 sup == with K(xt) -( 


min(x, t) J (11) 
1210,10) f, 91 G9dx 


max(x, t) 


Clearly, on I? [0, 1] the functional in (10) is defined everywhere except o = 0 (the zero 
function), and is bounded by 4/k, as proved before. Moreover, there is a clear homo- 
geneity property: The ratio for any q is equal to the ratio for any non-zero constant 
multiple cg, hence the ratio is constant on all rays {cp : c € R,c + 0}. 
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Therefore, the range of this quotient functional is clearly the same on the whole 
space L?[0, 1] V {0} as on B \ {0}, B denoting the unit ball B :- {p € L’[0,1] : lela < 1) 
(where the 2-norm ofa function in I? [0, 1] is lell, := (f, y’(x)dx)"?, as usual). Further- 
more, actually already on the unit sphere S :- (o € L2[0,1] : |l, = 1} the functional 
must take on all the values of its range. However, on the unit sphere the denominator 
is exactly one, so now we can modify the formulation and write 


11 11 
S(k) = sup | | g(x) o(t)K(x, t)dtdx = sup | | (x) p(t) K(x, t)dtdx. 
5050 B 0-0 


Moreover, it is clear that in this last formulation S(k) is taken by a maximizer function 
9 € L^[0,1] iff there is a maximum at some g € S iff there is a maximum on B (in 
which case again any maximum must belong to S). So any maximizer ¢ in the original 
formulation is maximizer together with all the ray (cq) of its homothetic copies, and 
in the new formulation this maximizer occurs exactly with c = +1/|lpl, i.e., at the unit 
norm elements of the given ray. 

This reformulation furnishes us the access to settle the existence question of some 
maximizer. In our formulation of the extremal problem all functions are real-valued, 
only for the next two propositions (spectral theory), and for the sake of being precise, 
we shall need complex-valued functions. 


Proposition 2. Let 
min(x, y) y 
max(x, y) 


K(x, y) := ( 
and define the Fredholm-type operator 


1 
T:L([0,1;C) > L’([0,1];C) (Tg)(x) := ILS y)dy. (12) 
0 


Then T is a compact, positive, self-adjoint operator on the complex Hilbert space 
L^([0, 1]; C), maps real-valued functions into real-valued ones, and preserves positivity. 


Proof. Since K € L®([0, 1] x [0, 1]), T is compact, see [9, 897]. Since 0 < K < land K 
is symmetric, the other two properties follow evidently. L1 


Proposition 3. S(k) is attained as a maximum by some maximizing function q € L’[0, 1]. 
Equivalently, 
-2 fo xq (x)q(x)dx 
fo xk-1q!2(x)dx 
for some appropriate q € X with q(1) = 0 and q # 0. 


S(k) = 


Proof. Consider the operator T as in (12). 


1: 3k 
Atv) = | | oWOOK(s ydxdy = Cro (13) 
00 
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which is a sesquilinear form on L?([0, 1]; C). By [9, 893] we have 
ITI = sup [Tol = sup (To, 9) - sup Ag, 9). 
Iell;si Iell;s1 lọl2<1 
Since T is compact, positive and self-adjoint, all of its eigenvalues are nonnegative, 
moreover, the eigenvalues can be ordered in a decreasing null-sequence (A), A, > 
ERA IE i A> 0 (j — co), and we also have 


IT] = max(A : Ais an eigenvalue of T} =: A,. 


Since T leaves the subspace of real-valued functions invariant, for any eigenvalue À € 
R of T there is a real-valued eigenfunction. Summing up, |T| = 4,, and there exists 
some (non-zero) eigenfunction q € L?[0, 1] satisfying ||g|| = 1 and A, = |T| = A(q, q), 
yielding a maximizer for A(g, q) as asserted. O 


Remark 4. The above proof yields also the following important information: S(k) is 
the largest eigenvalue A, of T, and any (normalized) eigenfunction ¢ of T belonging 
to A, is a maximizer; moreover, the only maximizers are non-zero eigenfunctions of T 
corresponding to A, = |[T |J. 

Indeed, as T is compact and self-adjoint, there is an orthonormal basis (ej) in 
L?[0, 1] that consists of eigenfunctions of T. Let 0 + 9 € I?[0, 1] be not an eigenfunc- 
tion to the eigenvalue A,. Then ọ = Yd (9, e)e; and 


oo 
(T9, 9) = X Ajo. e^ < Alel 
j=l 
by Parseval’s identity, where for the strict inequality “<” we have used that for some 
j > 1 we have |(9, ej] > 0, while 4; < A;. 

We also remark that since K is strictly positive, so is the operator T, hence one 
knows from Perron—Frobenius theory (see [7, Sec. 4.2]) that the dominant eigenvalue 
A, is simple with a corresponding strictly positive eigenfunction. This will be proved 
later also by directly determining all eigenfunctions of T. 


Next we turn to smoothness properties of eigenfunctions of T. 


3.2 Smoothness of maximizers and maximizers in C[0, 1] 


The above formulation also provides us a direct access to further smoothness state- 
ments. 


Lemma 5. The Fredholm-type operator T defined in (12) maps I? [0, 1] to the subspace 
C,(0, 1] of continuous functions with value 0 at 07. 


2 Notethat we identify functions defined on (0, 1] only but having limit O towards the boundary point 
0 with functions continuously extended to 0 by defining their value at 0 as 0. 
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Proof. Since L7[0,1] c L'[0,1], and 0 < K(x, y) € 1, the expression (Tg)(x) = 
IN 9(y)K(x, y)dy is an integral with a uniform majorant |g(y)| € L'[0, 1] of the in- 
tegrands. Hence by the Lebesgue Dominated Convergence Theorem, it suffices to take 
the pointwise limit under the integral sign. When x — x , this gives for all y > 0 
lm. K(x, vy) = K(x , y), while for y = 0 we have K(x,0) = K(x ,0) = 0 iden- 
tically. (Essentially, we have used only separate continuity of K on [0,1] x [0,1].) 
Thus lim, ,,,(Tp)(x) = f; p) limy- K(x, y)dy = f; (WK (xm y)dy = (Te) (x), i.e., 
Tọ € C[0, 1]. By definition, K(0, y) = 0 for all y € [0, 1], hence for every ọ € L”[0, 1] we 
have (T¢)(0) = 0. o 


Corollary 6. All eigenfunctions q of the Fredholm-type operator T defined in (12) are 
continuous and fulfill p(0) = 0. 


Equivalently, in the function space X defined in (7) all the functions q(x) = sf q (t)dt = 
- [in t ""o(t)dt corresponding to eigenfunctions q of T satisfy p(x) = x"?^q' (x) € Co(0, 1] 
and thus x"?-!a(x) e C(O, 1]. 


Proof. All eigenfunctions lie in the range of the operator T, hence belong to C,(0, 1] in 
view of Lemma 5. 


Recall that the correspondence between L?[0, 1] and our spaces Y and X was given by 
(x) = x"? p(x) = x"? q' (x). Thus for p € I?[0, 1], an eigenfunction of T, we obtain for 
the corresponding q that xl a (x) € C,(0, 1], whence also q € C(0,1] follows. More- 
over, lim, _,9, x”/?q'(x) = lim,_,, g(x) = 9(0) = 0, providing a continuous extension of 
x"!?q!(x) even to 0. Now writing q(x) = — [s q (t)dt yields q € C(0, 1]. While for x — 0+ 
we obtain that? 


1 
2: n[2-1 EE nf2-1f . -n[2 NE" 2 _ 
lim x q(x) = lim x ( [ox ar) = lim o(1) 2 0, 


x0 y= 


hence x"?-!q(x) € Cy(0, 1]. oO 


Although K is not everywhere continuous on [0, 1] x [0, 1], the operator T can still be 
approximated by compact operators given by continuous kernels. 


Proposition 7. The operator T restricted to C[0, 1] is a compact C[0, 1] — C[0, 1] oper- 
ator with exactly the same eigenvalues and eigenfunctions as on L?[0, 1]. 


Proof. Let f; : [0,1] x [0,1] — [0,1] be continuous with fixy) = 0 if x,y < z; and 
f(x, y) = Lif max(x, y) 2 " Then f;K is continuous, hence the integral operator T; 


n[2 


3 Hereo(1) means that for any e > 0 we have some ô such that for0 < t < ô, Iq (t) < et "^, Therefore, 


we have for any 0 « x ; ô the estimate |q(x) - q(0)] < E et "gt < b Um while x"?-! q(8) — 0 as 
n/2-1 2e 


q(x)| < <5 + o(1) and finally x" "2-1 a(x) > 0 with x > 04. 


x — 0+. Therefore, |x < 
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with kernel f;K is compact, see, e.g., [9, 890]. It is easy to see that T; — T in the 
operator norm (over C[0, 1]), hence T itself is compact, see [9, 876]. 


By Proposition 6 all eigenfunctions of T on I? [0,1] belong also to C,(0, 1], hence re- 
main eigenfunctions when T is considered as C[0, 1] — C[0, 1]. The converse is obvi- 
ous: Every continuous eigenfunction is of course also an eigenfunction from L7? [0, 1]. 
In particular, the set of eigenvalues are also exactly the same when considered in these 
two spaces. m 


Remark 8. One can show that the norm of T as an operator on C[0, 1] is 


n k: 


4 NR 4 NES 
ITlcii = (257) -(z) : (14) 


In fact, since T is positivity preserving, we have |[Tlcro;; = |T1/,,, where 1 is the con- 
stant 1 function. Easy calculation shows that (T1)(x) = 3 Qe xnl ? and that this 
function has maximum at x = (-27)/77, Since we already know |T|,2;0) = 41, and 
in general A, < ||Tllcjo1}, we obtain that the maximum of A(ọ, 9) with llọl = 1, i.e., 
S(K) is smaller than the constant in (14) above. 


3.3 Differentiability of maximizers 


We now push further the smoothness statements from the last subsection. We need 
some preparations, and define the following auxiliary functions 


mia) if(0,0) # (x, y) € [0 1]?, 
k(x, y) = max(x,y) ( ) (x y) | | (15) 
0 if (x, y) = (0,0), 
and Yuan ; n i 
2E i > > > 1 > 
T VEU if (0,0) + (a,b) € [0,1] (16) 
0 if (a, b) = (0,0). 
With these notations we have for every 0 < x, y < 1 the formula 
“K(x, y) = eG, y), K(x, y))kGo y); (17) 


which also holds for x or y being 0, with both sides vanishing. 

Now note that 0 < w(a,b) < nmax(a, p)" 2-1 < n. Furthermore, observe that for 
0 < a,b < 1 (even if both are zero) we have D"? — a"? = (Vb — Va) Xio aP = 
(b — a)w(a, b). Hence for any y, x, x’ € [0,1] we can write 


K(x', y) - K(x, y) = e(x(x', y), k(x, y)) - (K(x, y) - K(x, y)). (18) 
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Fix x » 0. Denoting A :- x' — x > 0, we also have 


<4 =-2 0 if y<x<x'; 
x x XX 
2 1 
1 E QUID a ; 1 
5 — > e ma etm m a f < < > 1 
K(x, y) — K(x, y) y F Ja uxsysx (19) 
E ue A if x « x' < y, 
y y y 
whence 
y y yA A A ; 
EE ee if y<x<x', 
2 
- A A : 
k(x’, y) - x(x, y)| = [-2|- 9227.22 ifx<ys<x', (Q0) 
A ae eee A : ] 
qz cy su if x«x «y. 
y 


Lemma 9. The operator T maps LI? (0, 1] to the space C'(0,1]. Furthermore, for p € 
D?[0, 1] and x € (0, 1] 


1 


1 
(Tg) (x) = J ez ke. y)dy = -(To)(x) + nx 7 | are (21) 
0 


x 


Proof. The two expressions given for (T9) (x) in (21) are easily seen to be equal, so the 
proof hinges upon showing that (TQ) (x) equals any one of them. 


For any 0 < x « x’ « 1 using (18) we can write 


1 
(Tg)x) - (Tox) _ ' K(x’, y) — K(x, y) 
aniline = | ey) w(k(x > y K(x, DD) dy. 


0 


(22) 


We fix x) > 0 and take either x = x, and x’ — x +, or x' = x, and x — x,-. In any 
case, by (20) and 0 < w < n we have the Lebesgue integrable majorant n|g( y)|/x, of 
the integrand, thus limit and integral can be interchanged. 
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For example in the case x x + taking into account (17) we are led to 


lim (Te)(x') - (Te)G) 


x! xg x! — xg 


1 
= | poale teo y) s 9) rne hoes 


0 


palrlo 3) KC» V) sgnQr— x) Seta 


Il 
o — 


1 
n 
Em | p(y) K(Xo; y) sgn(y — xo)dy 
Xo ; 


Xo 


1 
E AN o(y)K (xə y)dy - | P)KG y 


0 


1 
1 
= P(y)K (Xo y)dy - 5T Deo} (23) 


where we have used 2-k(x, y) = 2(x/y) = = 1/yfory > x and 2 s(x y) = x 2 (y[x) = 
-y/x for x < y. When: substituting the definition of K in the akow: we obtain all the 
asserted formulas. Note that in case y = xo, one-sided derivatives of x(xọ, y) still exist 
(and are equal to the limits from the respective side) but the existence and value of the 
limit at one exceptional point hes not interfere the value of the integral, therefore we 
have just put 0 for the value of 2 K(x, Ct hae , here. 

When x > x- = = x'-, the alstoni is entirely the same. 

The integrals on the right-hand side of (21) are of course integrals of integrable 
functions, and as such, are continuous in function of the limits of integration. There- 
fore, continuity of (T¢)' on (0, 1] also follows. oO 


Remark 10. When x, = 0, only the right-hand side derivative can be considered and 
thus we take x) = x = 0 and x' — 0+. Also, (Tọ)(0) = 0 and K(0, y) = 0, hence the 
consideration of the differential reduces to 


(Tp)(x!) 


1 1 gen 
e : Oe 
2 lim = | y(y)K(x', y)dy = Jim Lr D 
0 


lim 


x!'—0- x 


dy, (24) 


x! 


which, however, cannot be handled for general q € D [0, 1] or not even for g € C,(0, 1], 
and can be well estimated only if we use something more on q. See Corollary 12 below. 


Proposition 11. The operator T maps L?[0,1] to the space of absolutely continuous 
functions with bounded total variation. Moreover, for the total variation of Tp we have 


Proof. We already know that To € C!(0,1], so the total variation, whether finite or 
infinite, can be computed as V(T@, [0, 1]) = hs I(TQ)'|. Now the first formula from (21) 
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furnishes 
1 


1 
V(Tọ, [0, 1]) = | KT) (ldx = | 


0 0 


dx 


1 
ð 
| PO) ČK y)dy 


11 1 1 
ð ð 
< | [ione nara = [ieon( [e »|as )a < tens 
00 0 0 


as forall0 < y < 1 fixed we have IN |ZK(x, y)| dx = V(K(-,y),[0,1]) =2-y"? «2. O 


Corollary 12. If y lies in the range of T, then Ty € C'[0, 1], and (Ty)'(0) = 0. In partic- 
ular if ọ is an eigenfunction of T, then q € C! [0, 1], and g'(0) = 0. 


Proof. We have to calculate the limit in (24) for y := Tg in place of ¢. Recall from the 
above that then y = Tg € C,(0, 1], in particular y(0) = (Tg)(0) = 0, and y € C! (0, 1], 
V(y,[0,1]) x 2lgl,;. The second mean value theorem and integration by parts yield 
with some appropriate z := z, € (0, x) 


1 : I m[2 
L | W(y)K(x', y)dy = | | w(y)K(x's y)dy + | d 
0 * 0 s 


x! 
x! 


In[2-1 /2- 


-x 
n[2-1 


= w(z,)K(x,z,) + {| vo 


1 L x" i 
1-n/2 ! 1-n/23 
y ML E l, 


so the first term tends to y(0) = 0 when x’ — 0+. The term in the square bracket 
m[2-1-yQ) yx) 
n[2-1 n[2-1? 


converge to 0. Finally, for the integral Proposition 11 gives that y’ € L'[0,1], while 
the product of the further factors stays bounded uniformly for all x’, y € [0,1], as the 
integral runs only through values y > x’. That is, we can again use the Lebesgue Dom- 
inated Convergence Theorem and calculate the limit by moving it under the integral 
sign. Furthermore, the pointwise limit of the expression is zero for all fixed y, whence 
the assertion follows. o 


contributes x 


and as x’ — 0+ and w(x’) — y(0) = 0, both terms 


4 Solving the maximization problem 


4.1 Setting up a differential equation for potential extremal 
functions 


By the previous section we know that our maximization problem has a solution, and 
we also saw that maximizers are sufficiently smooth. We can now set up a differential 
equation to find maximizers, or which is essentially equivalent, to find the eigenfunc- 
tions of T. 
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Proposition 13. Let o € L”[0, 1] be an eigenfunction of T corresponding to the eigen- 
value à > 0. Then g is continuous on [0, 1], infinitely often differentiable on (0, 1]. The 
function q(x) = - lk gly) yr dy satisfies q(1) = 0 and the differential equation 


n n b E es n 
q'(x) + <4 (x) + =4(x) -0 (xe(01]), where b:= i^ 0. (25) 


Conversely, let A > 0 and suppose that q is a non-zero, C^ (0, 1] solution of the differen- 
tial equation above with q(1) = 0. If p(x) = xl Ed (x) extends continuously to 0 with 
lim, ,9, x”/q'(x) = 0, then q is an eigenfunction of T corresponding to the eigenvalue 


À» 0. 


Proof. If ọ € L?[0,1] is an eigenfunction of T for the eigenvalue A > 0, then it be- 
longs to the range of T, hence is continuous and continuously differentiable on (0, 1] 
by Lemma 9. Substituting Tq = Aq in (21) we obtain 


wae 2 Aun pt [22 9) 5 d 


n[2 
xo 


As the right-hand side is differentiable, we can differentiate also the left-hand side 
showing q € C^(0, 1]. We substitute x”? p(x) = q'(x) and g(x) = x"?^q' (x) and obtain 


1 
d i Nyon n[2- 
3; 0x74 09) = sag (xe) + nx [d oa, 


and hence j 
zm (Ax""q'()) =- L AzX x"? 1 q'(x) - nx"! a(x). 


Differentiation yields 
Ax"? qo) + AZ! P7 qr (x) = - T aq (9) ~ nx"? g(x), 
and then by rearranging we obtain 
Ax"? q" (x) + àz "e? q (x) nx"! g(x) =0 


n/2 


Division by Ax""^ thus leads to the asserted differential equation (25). 


To see the converse we set y = Tg. Note that then y € C,(0, 1] according to Lemma 5. 
Then q(x) = x”/?q'(x) entails 


-n[2-1 


q'(x) = ES g(x) + x "^g 


x), 


so that using the assumption that q solves (25) we obtain 


1 
N -n/2- -n -n n/À -n 
zs m *p(x) px 9 (x) + zx (x) + nih -= IE gody) =0, 
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and thus also 
1 


TAr" gla) + Ax" g'a) - z f y eondy = 0. 


By (21) with y = Tg we also have 


1 
x"? I x "2 1 n AY) y 
y'(x) = voi y, 
so 
Sx 9(x) 4 Axe n/2 9 (x) - x"? y (x) - zx yx) 20 


If we multiply by x”, we obtain for all x > 0 
0- z Ax gla) + Ax"? gl (x) — x" Py! (x) — x" w(x) = i (x" Ap(x) - x”? w(x). 


Since p, y € C[0, 1], x"? (Ag(x) - v(x)) must vanish at 0, whence Ag = y follows. O 


Thus the solution of the maximization problem is reduced to solving the homo- 
geneous second order ordinary differential equation (25), and to finding the feasible 
values of A. Next we solve this equation and analyze some properties of the solutions. 


4.2 Bessel functions and Bessel’s differential equation 


Recall Bessel’s differential equation 
[D 1, y? 
y (x) + cy GO c 1-5 jy) 70, (26) 


for some fixed parameter v € IR. The Bessel function J, (of the first kind) is a solution 
of this equation, [2, (6.71), p. 115]. Notice that J_, is also a solution, but for v integer 
J, and J_, are linearly dependent, in fact J_, = (-1)”J,. So for v € N another, linearly 
independent solution is needed, which is provided by the Bessel functions Y, (of the 
second kind), obtained as 


Y,(x) = he) | a, 
where specifying the limits of integration is equivalent to fix some primitive of the 
integrand, and the integration limit cannot be at 0 where Y, (x) is divergent in the order 
x”, see [2, (6.73), (6.74)]. Then for v € IN the general solution of equation (26) is a linear 
combination cJ, + c; Y,, see [2, 8102]. 
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4.3 Solution of the differential equation 


Bowman computes, see [2, (6.80), p. 117], what happens if we consider the trans- 
formed, substituted functions u(x) := x" y(fix"), where y satisfies the Bessel equation 
(26), and establishes that then the new functions u(x) will be the general solutions of 
the transformed equation 


u! (x) - 2a 


2.22 
Ly (x) + (ev aye MJ u(x) = 0. (27) 


If we choose here the parameters v := m, a := -m/2 = 1-k/2, B := 24n/A and y := 1/2 
(where n :- k - 1 = m+ 1 and m := k - 2 as fixed above in (8)), then the equation 
(27) becomes exactly (25). Thus we obtain that for any fixed values of m := k - 2 and 
À » 0, there is a one-to-one correspondence between the solutions q of (25) and y of 


(26) given by q(x) := x“ y(fix") = x! ^ y(2V/n/Ayx). 
Corollary 14. Every solution q of (25) is a linear combination of transformed Bessel func- 


tions from the above, i.e., 


q(x) = qx "7, (2VnJAVx) ex ""Y, QNnlA Vx) (m=n-1=k-2). 


4.4 Some analysis of the occurring Bessel type functions 


Before proceeding, let us note one important thing. Not all solutions of the differential 
equation are relevant for us, because the resulting q’ must have finite weighted square 
integral, i.e., p(x) = xg q (x) has to belong to I?[0, 1] (and even be continuous on 
[0, 1], according to dicus 13). 

According to the second formula of [1, 9.1.30] (with the choice k = 1 and v = m 
there) we have 


ld = 
(4) (x "aJ, (X) + GY, ,(x)) 
x ae "(erasa (x) + ©Yinei(X)) = E^ (eJ, G) + oY,(x)), 


therefore we obtain for a solution q of (25) that 


q'(x) = Z (ax om? QNIN) + x? Y,,(2Vn/AVx)) 


i = x "P (op (NAAND + oY (a RIA). 


So that 


n/2 q( 


g(x) = x""^q'(x) = (cJn(2Vn/A-VX) + QY,(2-Vn/A Vx) . (28) 


Vx 
ET: 
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To see when such a ọ may actually be an eigenfunction of T, we first find out when it 
belongs to C,(0, 1]. First, see [2, (1.2)] 


e CD x 


— v — ~ v 
Jt x 2 Tu AQ whencealso J,(x) ~ =x" (x > 04), (Q9) 
so in particular J, is continuous on [0, co), and for any e, € R the part LEG Jn 


QA n/ A Vx) belongs to C,(0, 1] c L’[0, 1]. Second, Y,,(x) = x^" (x — 0+) [2, p. 116], en- 
tailing that forc, #0 oY, (av) = x")? and this function is not even bounded 
near 0, ifn > 1 (i.e., when k = n-1 > 2). Ifn = 1, i.e., k = 2, this function is not 
vanishing at 0, a condition that is necessary for an eigenfunction of T by Lemma 5. 
So from Corollary 6 we obtain that q in (28) belongs to C,(0, 1] (and thus may be a 
candidate for being an eigenfunction of T) if and only if c, = 0. 


As a consequence of this and of Proposition 13 we obtain the following. 


Corollary 15. Consider the operator T from (12), and let à > 0. Then à > 0 is an eigen- 
value of T if and only if 
Im (2Vn/A) = 0. 
In this case 
g(x) = x" ^g (x) = cx ^T (2VnjAvx), qa #0 


are the only eigenfunctions corresponding to À. 


For given m let us denote the roots of J,, by a,,,.(r € IN) ordered increasingly. At this 
stage it is in order to recall the following about the zeros of Bessel functions. We have 
Amr — œ (r — oo), and for rather large values the roots œ, of J,, are very well 
distributed, as essentially there falls one root in each interval of length z. However, 
for fixed m the increasing sequence of zeros (æ+) starts only with æn, ~ m + cm! ^, 


with c = 1.8557571 ..., see [1, 9.5.14, p. 341]. Let us introduce the notation 
EP EEECCERVTM 


and 


el) = qu) = x "PT 29] M0 + Dany x) = x "P T, (Cn X). 


By putting everything together we obtain the following result. 
Theorem 16. For the extremal problem (1) we have 
| 4(k- 1) 4 


S(k) =A, = —— = ———— 
TUM ay, — k*2ck!5 +O(1) 


(30) 


with o,_,, the first root of the order k — 2 Bessel function J,_,(x), and the constant c = 
1.8557571.... 
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The only extremal functions for the formulation (5) are non-zero constant multiples 


of 
qi) = gcnus p (Apa vx). 


Proof. Uniqueness follows from Remark 4 and Corollary 15. For q, being a maximizer 
for (5) it remains only to show that 


1 
[aora +0. 
0 


But this follows since q; is (strictly) positive all over (0,1), a,_,, being the very first 
zero of the Bessel function J; ;. o 


5 Power series 


The above settles the issue ofthe best weight P — and also the order of S(k) — in the GPY 

method. However, not all weight functions are easy to handle, and a Bessel function — 

even if an analytic function with relatively strongly convergent series expansion — may 

beunmanageable, at least in our current technical abilities. We will discuss, using the 

classical series expansion of J,,, how it may work in this context. Actually, not too well. 
With the notations from the above for q,, o,, , etc., (29) yields 


2 


| .-m[2 g-m ES (-b)/ j ; A _ g2_ mr 
q(x) = x "^g WORNI Tam with B:= Vb, b= = T 


Then we can evaluate the functionals 
1 1 1 
G(q,) = | x""g"(x)x, F(q,) = fon + 1)x"q'(x)dx = (-2) | x"*!g (x)qi(x)dx. 
0 0 0 


Computations with the power series provide 


"DEDE" M 
T o fe mitte Djs je Dn t DU 
and 
co (-by** 
A MR M 
F(q,) EA (m4 j«£«2)lém s Dm e+ 1) 
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With some reformulations we can also write 


co (-b)” v 1 
F(q) 2 2:56 mv) (m+v+2)! È j= (CETT ) 


G(q,) b xe a CU a y% mE) 


v=0 (m+v+1)!(m+v+2)! &j=0 \ v-j 
oo (=b)” 2m*2v4l co . (-b'Qme-vi2) (2m*2v«2 
x 2 Mio aa v ) - 1 25-0 ERENCE v ) 
= oo (Cb Im+2vI2) bye (Cb 2m+2v+2 
b 25/50 omm v ) b cd ——— v ) 
co | (-bD' (2m+2v+2\(2m+2v+1 co |  (Q2ms2vel) —— O pv 
i 2 Cmca m+v X v ) aA dvb (EX amil b) 
~ by (D (2m+2v43)(2m42v42) b5 —— Omt?! py 
208 Were m+v+1 X v ) 2550 Gaia DOMED b) 


Unfortunately the series expansions here have large and oscillating terms, so deal- 
ing with it does not seem to be simple. When, e.g., b is of the order m?, then also the 
terms with v « m are the highest, and there are a large number of similar order large 
terms. Therefore, this series expansion does not seem to be suitable neither for the 
computation of the value of the ratio, nor for the extraction of a good polynomial ap- 
proximation which would approach the global maximum while remaining manage- 
able. 


6 Approximate maximization by polynomials 


In the aimed applications in showing small gaps between consecutive prime numbers 
it is very important to have a suitably nice, manageable function P. It is enough to 
mention that even in the simplest case of P(x) = xE, £ = vk the technical difficul- 
ties become rather serious when k and £ tend to infinity with the size N of the primes, 
see [5]. The details of these aspects, when the choice of the weight function is done 
according to the present work, will be handled in the forthcoming paper [8]. Here we 
will only present the foreseen choice of the weight P, and show its approximate opti- 
mality. The said choice will be a relatively simple function, actually a real polynomial 
P(x), satisfying the conditions 


x"|P(x), P(1) » 6 deg P(x) - k+ Cok", 31) 


which is essentially optimal in the extremal problem (1). More exactly, with the nota- 
tions 


l ka l kl 
Ay = | a— (P*Y1-x)) dx, B, =|= " (P?a-x))*dx (2) 


0 0 


it satisfies with an absolute constant C, 


A,(k + C,k'/?) — 4B, > 0. (33) 
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Equivalently, 
Ak 4 


By im C, ki" 
in full correspondence with (30). Compared with the optimal transformed Bessel func- 
tion q,(x) = x nn Jin (%m,1 V X), the difference is only in the value of the constant C}. 

In order to define our polynomial we put 


M :=[C,k'7/6], g(y):- (y- D*2- yy’. (35) 


S(P, k) := (34) 


Let us remark that the exact choice of g( y) is irrelevant, any positive polynomial or 
even a function g € C'[1,2] with a zero of order at least 3 at y = 1 and y = 2 would 
suffice for our purposes. After this, let 


P(x) :- P(x) = Yo (£) (5y on . (36) 
m 


In evaluating A, and B, we will use the well-known relation (easily obtained by partial 
integration and induction) for the Euler integral 
1 1 
B(m,n) := | x"(1- x)"dx = MN... NR | x dx eee (37) 


(n+1)---(m+n) "(nen Dd! 
0 0 


In view of (36)- (37) we have 


PN Ye«(£) (5) i (39) 
241 285 
eney- Y YX)G0() ze a 


A= Y Y a(S) o(2) BC D) (42) 


RUNE E LA (kI pe e£ 
sal uui e k 
k 22.8 M/)9\M/) (k*& t & (3) 
2H] 2469 
In the following, put 
u:-£ 1 v:6, +1, and H := (u+ v{v)/2, D := (u-v)/2. (44) 
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Taking into account 2 + £,, £, and g(1) = g(2) = 0, the even variables u, v will run from 
M + 1to 2M and we will have 


Ketek M? 
— = 1+0| — |. 45 
(k+£ +6)! n ( k ) m 


Here and elsewhere in the sequel the implied absolute constants of the O symbol as 
well as the absolute constants C; (i > 2) will be always independent from C}. 

Clearly, we can replace A, := k!A, and B; = k!By for A, and B,, resp., in (33), 
hence in order to show (33) it suffices to prove 


2M M? 
PEU) 20, (46) 


where 


E u-1 v-1\/u+v0-2 k * CK? (u+v)(u+v-1) 
(H) = 2" » o M Ja M X u-l1 ) EXAM uv E 


2|u,v, u+v=2H 
u,ve(M,2M] 

(47) 

Let us denote the above summation conditions simply by Y, ur and let us consider first 


L(H):- py E wie) [ERE al 
s25" a (#1) (“te 7) 2a TE 


We will see later that the mere significance of the weight function g is to cut the tails 
when H is near to M or 2M, so first we will investigate the simpler sum 


n Hy (u*v-2| 4D! -2H 
gane ry (pn “2 
when 
H € [M +t(M),2M -t(M)], t(M) := 4Mlog M. (50) 
By Stirling’s formula 
log I'(s) = (s - i) logs - s + $ log(27) + o(4), (51) 
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we obtain 
3gíu*v-2 I(2H -1)2?H 
eee?) aye apenas: 
u-1 I(H + D)I'(H - D) 


log 2 
(an - :) (log H + 1og 2 + log (1 - =) +1- 98 28 - 2Hlog2+0(=) 
2 2H 2 H 


- (H+ D~ +) (logH +1og(1+ 2))- (H-D~ 2) (log H + tog(1- 2)) 


_ riogH _3log2 || , _ logan 
2 2 
co 2D D 2ntl 2H D 2n+2 1 D? 
- Mi) eee )+0(5+2) 
_-logH logl6oxr <= 1 pres 1 D 
RE TA -2 Daam (zm) 


Using e! = 1 + O(x) we obtain from (49), (50) and (52) 


1 E D? I D D 
e a Y £9" (reo z mm 
2JA4H9? A, H H'H'H 


M<H-D,H+D<2M 
[o9] 
(H 1 Ü _ I; (H 1 
= PUR ) «o()«o( | —e Ua) - EG ) «o(). 
2 VAH?’ M? H IVAH?’ M? 
tM) 


where 


co 2 j 
I, (H) := 25 f(m), f(x) = e = je" IH 


m--oo 


We will use Poisson summation formula 


> f(t + nT) = 7» AZ) e2mivtlT 


n=—OO 


with T = 2, t = 0 for 2|H and T = 2, t = 1 for 2 + H, valid if* 


Lf G9] + Ifl <C(1+ lxp[ 173) with some à » 0,C » 0. 


4 For this form see [12, Chapter VII, 82]. 
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(54) 


(55) 


(56) 
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Let us consider first the case when H is even, that is when m - 2n in the above 
sum (54). We have then 


2I (H) -2 Y. fQm)- Y f(vi2- Y, | fle" dx 


n-—oo v=-00 v--oo | 
co co 
V -x?/H x 7ivx S -y 2 nivVH 
=> IE 2— -1]e dx- NH Y [eG -9. ?dy 
v=- $5 H v--oo _ 
E VH X g Uri HY m v H/4 
v--oo 


. {ay - nivVH/2)° + 2(y- nivVH/2)nivVH — m H/2 - 1} dy 


oo 
= VH > | ee 127: + 2nivNHz - m Y! H[2 - 1} dz 
TOPs SERVER 


oo 


= VH » | QU HH bz + 2nivNHz - n y H[2 — 1} dz 
VECO 3220 
= VH | e” (22 - Ddz + O (e?) - VH [-*]^ + o (e°) = o (e°). 
—oo 
zeR 


(57) 


The proof runs completely analogously fort = 1, that is, when H is odd, because the 
extra factor e™ does not change the modulus of f (v/2). 

Dealing with the original integral I; (H) we have to take into account the effect of 
the weight function g as well. From the Taylor expansion of g we find 


HxD-1V (H-1 i(H-1\D p 
(Uu) a cm 
so this effect is 
H+D-1\ ÍH-D-1V . ;(H-1 D’ 
a a E — 09 


The effect of the error term on (H ) is here, similarly to (49), 


2 4 2 

o(r Y ges + al + a: ) =) 5 (=) re 
2|D-H A H M M 
|D|<M 


So we obtain for all integers H subject to (50) from (48), (49), (53), (55), (57), (59) and 


(60) 
a( ma 
=) e(5) C5 


L(H)- wang (53) * oss = Tepe 0D +o(st) 2-55 (61) 
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On the other hand, if (50) does not hold, that is if 
min (H - M,2M - H) < t(M), (62) 


then we find [D| < t(M) and |H € D- 1| € [M, M + 2t(M)]U [2M -2t(M), 2M], whence 


H+D-1 H-D-1 t*(M) 1 
( M ) ( M )< MS ME (62) 
This implies in the case of (62) 
foe} A ae 3 
(H) « M7 MP | e” E| Z +1) du«w ^. (64) 
k H 
Consequently, for k > ky we have by I, (H) > I,(H) 
)Ó nis X  Ļ®)2-0M”. (65) 
HOM MH sic Manet 
Finally, let us investigate now for H in (50) the difference 
M .— 
1,(H) - L(H) = q KD, (66) 


where similarly to the above considerations 
Pe p * -1 v-1\/u+v-2 
a sre) 
(H) Lad ar 9r | ui 


1 -D'iH ( ( 1 D? D’ ) 
= e 1+ O = SS o—— 67 
2JRH 2. H P Hm (67) 


M«H-D,H*«Dz2M 


Summing over all H in (62) we obtain 


2M-t(M) 2M-t(M) , C M? 
Y He Y H+ p (68) 
H=M+t(M) H=M+t(M) 


Therefore by (61) and (65)-(68) we have finally 


2M 2 
CM C C 
Ae Se >o, (69) 
H=M+1 k M M 
if C, was chosen sufficiently large, satisfying 
EVE 
C,>6 (2) . (70) 
Cs 
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7 Note added in proof on June 28, 2013 


This work was done during the academic year 2012/2013 with the aim of finding the 
limit of the GPY method (if the method is applied without any further improvement 
in the admissible distribution level of primes in arithmetic progressions, i.e. without 
any further improvements of any Bombieri-Vinogradov type theorem, as was the case 
a few weeks ago). However the situation changed during this short period of time. 
The recent breakthrough of Y. Zhang, Bounded Gaps Between Primes, Ann. of Math., 
to appear, reached a higher level of distribution of primes than the earlier value 1/2 
in smooth moduli for special reduced residue classes. This can substitute a uniform 
improvement of the Bombieri-Vinogradov theorem in this application, as shown in 
the work of Y. Motohashi and J. Pintz, A smoothed GPY sieve, Bull. London Math. Soc. 
40 (2008), no. 2, 298-310. In such a way Zhang could show that there are infinitely 
many prime-pairs with a distance of at most 70 000 000 apart. The application of the 
Bessel-kernel now helps to significantly reduce this bound. Using this and many other 
important ideas the Polymath project of T. Tao (see *The prime tuples conjecture, sieve 
theory, and the work of Goldston- Pintz-Yildirim, Motohashi-Pintz, and Zhang"... and 
other parts of Polymath 8 at wordpress . com) succeeded to diminish this bound below 
7 000. In this context we remark that the careful investigation of the problem (the ex- 
act proof of the existence of the maximum of the quantity S(P, k) of (34)) in a very 
wide class of functions, furthermore its substitution with an easier manageable poly- 
nomial in the last section of the work) are not necessary for the improvement of the 
bound. However, apart from the improvements in the size of the gap, these consider- 
ations show that there is no room for possible further improvements in diminishing 
Zhang's bound if one uses a different kernel function. According to this, the unpub- 
lished investigations of J. B. Conrey and his colleagues at AIM which first found the 
Bessel function in 2005 could have played the same role. 
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